Rashba Spin Orbit Interaction and Birefringent Electron Optics in Graphene 
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Electron optics exploits the analogies between rays in geometrical optics and electron trajectories, 
leading to interesting insights and potential applications. Graphene, with its two- dimensionality and 
photon-like behavior of its charge carriers, is the perfect candidate for the exploitation of electron 
optics. We show that a circular gate- or doping-controlled region in the presence of Rashba spin- 
orbit interaction in graphene may indeed behave as a Veselago electronic lens but with two different 
indices of refraction. We demonstrate that this birefringence results in complex caustics patterns 
for a circular gate, selective focusing of different spins, and the possible direct measurement of the 
Rashba coupling strength in scanning probe experiments. 
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The analogies between geometrical optics and electron 
trajectories have resulted in a number of interesting pro- 
posals for device applications pQ, where interfaces play 
a similar role to that played by transparent interfaces in 
physical optics. This leads to the manipulation and con- 
trol of electron trajectories, where the major factor de- 
termining the electron dynamics is the change in group 
velocity through these interfaces, mimicking refringent 
physical optics and lenses. Such change in group veloc- 
ity is typically achieved by local gating, which modulates 
carrier densities and fixes the corresponding index of re- 
fraction. 

Optical birefringence in materials results from crystal 
anisotropics which are manifested as different group ve- 
locities for different polarizations of the propagating light 
in the material. In this paper, we show that an equivalent 
phenomenon to optical birefringence in electron optics is 
feasible in two dimensional graphene, which in essence 
reflects the intrinsic crystal structure even at large elec- 
tronic wavelengths. The effect requires the presence of 
Rashba spin-orbit interaction, where the different group 
velocities depend on the chirality of the electronic states, 
mimicking the light polarization dependence of the group 
velocities in optical birefringent materials. 

The low energy dispersion of electrons in graphene is 
centered near two inequivalent points in the Brillouin 
zone, the K and K' or Dirac points [2H4]. The "mass- 
less" nature of electrons results in novel phenomena such 
as the Klein paradox [3j [5j [6] , which leads to full trans- 
parency of a sharp gated interface for normal incident 
electrons, and a high probability of transmission for in- 
coming electrons with finite angles. The linear disper- 
sion of electrons is also evocative of photons, prompting 
a number of proposals and experiments to probe optical 
analogs with charge carriers [7], aided in great measure 
by the high electron mobility in this unique material. In 
fact, the transparency of barriers and ability to gate re- 
gions of the system to change the sign of carriers can lead 
to the use of graphene gate-controlled interfaces as elec- 
tronic lenses that follow Snell's law with negative index 



of refraction and allow the implementation of electronic 
analogues of Veselago optics [8HTD] . 

The study of spin transport properties of suspended 
and deposited graphene is rapidly becoming an impor- 
tant area of research. Experiments have achieved spin 
polarized injection of electrons and measure spin valve ef- 
fects [TTHT3] , spin polarized currents with long coherence 
lengths in suspended graphene [H] . and shorter coher- 
ence lengths for deposited samples or samples containing 
impurities that enhance spin orbit effects [HI [16] , such 
as hydrogen or gold [ill EH] • An important ingredient 
determining the carrier spin dynamics in this and other 
materials is the spin orbit interaction (SOI). The intrinsic 
SOI respects all the lattice symmetries in graphene and 
results in a small energy gap at the Dirac points [19] . The 
extrinsic or Rashba SOI results from the lack of inversion 
symmetry due to perpendicular electric fields, substrate 
effects, chemical doping, or curvature of graphene cor- 
rugations [20j [21] . In this paper we focus on the effects 
of the Rashba SOI in the spin dynamics of carriers away 
from the neutrality point as the spin-orbit interaction can 
be controlled by external factors [22] . 

To study these effects, we consider a graphene sheet 
with a circular gate potential (or corresponding dop- 
ing/intercalation profile) covering an area of radius 
where the electric field (or doping) reverses the car- 
rier character from electron to hole and also generates 
a Rashba SOI of strength as shown in Figure [lji. 
Electron scattering from such gated regions has been 
shown to result in pronounced wave function intensity 
maxima due the constructive interference between dif- 
ferent partial wave components [8]. These patterns are 
similar to the optical caustics which develop by light re- 
fraction through a shaped medium, and belong to a class 
of cusps in catastrophe theory [23 . We should comment 
that recent experiments have demonstrated gate tunabil- 
ity, which makes the study of electron optics on graphene 
viable and controllable [6j [24] [25] , Moreover, imaging 
of the patterns resulting from the electron flow through 
these nanoscale structures can be achieved through STM 
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FIG. 1. (Color online) a) Electron flux in graphene coming along the x-direction onto a circular gate potential covering an area 
of radius R; the normal electric field reverses the carrier character from electron to hole and also generates a Rashba spin-orbit 
field, b) Three dimensional map of \i/j\ 2 near and inside gated region centered at (x,y) — (0,0) (normalized to incident flux). 
The electronic patterns produced inside the gated region are analogous to those produced by the refraction of light in a circular 
medium with refractive index n = — 1 (Veselago lens). The electronic index of refraction is determined by the ratio of wave 
numbers in and out of the gated obstacle: both regions are characterized by a linear electronic dispersion (upper left inset), 
and in this case, kR = 300 (electron) and k' R — 300 (hole). High intensity maxima correspond to caustics with p = 1 (x < 0) 
and p = 2 (x > 0), where p — 1 is the number of internal reflections, x and y coordinates in units of the radius R. c) In the 
presence of Rashba spin orbit interaction the electronic dispersion in the gated region is modified, so that scattering particles 
have access to two different wave numbers (upper right inset). Scattering produces electronic patterns as shown in the three 
dimensional map \ip\ 2 , with kR = 300, VR/Hvf = 600, XrR/Hvf = 100, for "[-spin incoming flux. The p = 1 and p — 2 caustics 
and cusps are doubled, associated with states of different chirality in graphene. The pattern can be described by the refraction 
from a medium with two (negative) indices, n±. Here, the degree of birefringence is An = n_ — n+ = 0.71. The dual pattern 
persists even for spin-unpolarized incidence. 



[26j [27] or other scanning probe techniques [28] [29] . 

The Hamiltonian of the system has the form H = H Q + 
Hy + Hr, where H /hvF = cr x k x + a yky describes free 
electrons in graphene with momentum k — (k Xl k y ) away 
from the K point. Hy = V$(R — r) represents the gated 
or doped region of strength V and Hr = \R(cr x Sy — 
s x a y ) , &(R — r) is the Rashba SOI in this circular region 
for k near the K point [4j [19], where {cr^} and {s^} 
are Pauli matrices representing the electron pseudospin 
(A,B) and spin (t,|), respectively, A# is the strength of 
the Rashba interaction, and $ is the Heaviside function 
determining the radius R of the region. The Hamiltonian 
is parameterized in terms of hvp, where vp = 10 6 m/s 
is the carrier Fermi velocity [2j [3]. We are interested 



in Rashba regions with R ^> a, where a is the lattice 
constant of graphene, in order to adopt the continuum 
description of graphene. 

The cylindrical symmetry of the system allows one 
to write the eigenfunctions in terms of four linear dif- 
ferential equations coupling the different spinor compo- 
nents, i/jj (^ A ^^ B ^^Ai^Bi) T The symmetries re- 
sult in the conservation of the z-component of the total 
angular momentum in the system, [J Z ,H] = 0, where 
J z = L z + hs z /2 + ha z /2, and L z = —iftde is the orbital 
angular momentum; as such, the eigenstates describing 
the system can be labeled according to their total angu- 
lar momentum, j, with J z ipj = jhipj (see Supplemental 
Material). 
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The scattering of an incident flux of electrons on the 
gated obstacle is studied through a spin-dependent gen- 
eralization of the partial wave component method [30 . 
This utilizes the analytical solutions for the spinors to- 
gether with their asymptotic expansions and bound- 
ary conditions in order to obtain the scattering ampli- 
tudes. The latter allow the direct evaluation of wave 
functions inside and outside the Rashba region, depend- 
ing on the incident energy, momentum and spin con- 
tent of the incident flux, as well as system parameters 
(see Supplemental Material). The momentum outside 
the scattering region is given by k = E/hvp, while 
the momentum inside takes two different values, k± = 

a / (E - V) 2 ± 2\ R (E - V) /hv F , for the given energy E, 
due to the presence of the Rashba SOI (see Figure [T]). 

Formation of caustics has been described in the ab- 
sence of SOI [10 , resulting in beautiful scattering pat- 
terns. The caustics and cusps can be understood in terms 
of geometrical optics, and characterized by a negative in- 
dex of refraction n = —\E — V\/\E\ = —k'/k, where k' 
is the wave number inside the gated region [8 (Vese- 
lago lensing [9]). For an incident electron flux along 
the x-direction, the position of the cusps produced in- 
side a circular gated region can be shown to be given by 
%cusp — (— l) p /(l n l — 1 + for p — 1 internal reflections 
of the ray inside the region [10] (see Figure [TJd). 

In contrast, in a system with Rashba SOI the two wave 
numbers for a given energy allow for two refraction in- 
dices, associated with the two chiral solutions of the Dirac 
equation, n± = —k±/k. Correspondingly, the optical 
analogue results in two different cusp locations (see Fig. 
[I};) given by 

*± = ( " 1)P (1) 
cusp |n±|-l + 2p ' 1 ; 

As we will see below, this simple description is borne 
out by the full quantum calculation of the scattering pat- 
terns. We should emphasize that despite the optical anal- 
ogy, birefringence here has a clear quantum mechanical 
origin, as the group velocity and built in phases of the 
different chiral states are different due to the presence of 
the Rashba interaction. 

The presence of the Rashba interaction causes the spin 
of the incoming electron to precess as it travels along 
the scattering region; this spin precession is seen as os- 
cillations in the amplitudes of the spin components of 
the wave function (Rashba oscillations), where a clear 
spatial frequency can be identified for different A# val- 
ues: shorter precession wavelength for larger A#. For 
small Rashba coupling, the wavelength of these oscilla- 
tions is comparable to the size of the gated region, lead- 
ing to the presence of different spins over large areas of 
the scattering disk. Figure [2^i,b shows the wave func- 
tions inside and outside the scattering region for spin 
up and down components for \rR/Tivf = 3; take no- 




FIG. 2. (Color online) Probability density patterns (scale bars 
show log of amplitudes, normalized to incident flux) resulting 
from the scattering of incoming \-spin electron wave along 
the x-direction with kR — 300, in the presence of a gate 
potential VR/Tivf — 600, and Rashba coupling XrR/Kvf = 
3. This arrangement leads to a system with a small degree of 
birefringence An = n_ — n+ = 0.02. a) The spin-preserving 
component \ip^\ 2 an d b) spin- flip component \^i\ 2 display 
Rashba oscillations with a wavelength covering a large region 
of the scattering target, due to the small A^. c) Total wave 
function \ip\ 2 = \ip^\ 2 + l^-ul^ displays caustics and cusps 
that remain almost unchanged from the case in which the 
Rashba interaction is absent, d) Net spin n ~ l^ttl 2 — iV'til 2 
quantifies the predominance of a given spin in different regions 
of the system. The large wavelength of Rashba oscillations 
leads to a cusp at x ~ —0.5 with a net ^-spin (blue) while the 
annular region near x c± has a net t-spin (red). 

tice these are results for spin-up incidence. Notice that 
both spin components display caustics inside the scatter- 
ing region. Moreover, one observes that the wave func- 
tion around the caustics has a net spin content, as shown 
by rj = sgn(A)|log 10 (|A|)|, where A = |^ n | 2 - l^ul 2 > 
in Figure [2]i. For this set of parameters, the p = 1 
caustic (x < 0) is predominantly spin-down, while the 
p = 2 caustic (x > 0) is predominantly spin-up. The 
total wave function, \i(j\ 2 = \ip^\ 2 + |^nl 2 snows no ma- 
jor changes due to the Rashba SOI, since A# is small 
and the two different wave numbers are nearly equal, 
k- — fc + ~ 2\r/Tivf] see Figure [2^. Correspondingly, 
the two indices are also very similar, n_ — n + = 0.02, 
and the pattern of caustics and cusps is essentially un- 
changed from the A# = case. 



For larger values of the Rashba interaction the oscilla- 
tions become even sharper, with shorter wavelength, and 
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FIG. 3. (Color online) Wave function maps (scale bar show log 
of amplitudes; normalization to incident flux) similar to those 
in Figure[2]but for larger X R . Here, kR = 250, VR/Tiv F = 600 
and XrR/Hvf = 30. This arrangement leads to a system with 
a degree of birefringence An = n_ — n+ = 0.24. The different 
energy of the incoming wave leads to modified positions of 
caustics and cusps, a) l^ttl 2 an d b) iV^vlJ 2 snow clearer and 
sharper Rashba oscillations with shorter wavelength. Notice 
that the number of oscillations inside the gated region is ap- 
proximately XrR/2Tivf — 15 for both spins, as shown by the 
rings around the main cusp, c) The total wave function, l^l 2 , 
displays clear duplicate sets of caustics and cusps manifesting 
the birefringent character of the scattering, due to the large 
Rashba interaction and associated different group velocities, 
d) The net spin 77, shows variation of the two spin components 
along the the different cusps and caustics. 



their number inside the gated region is approximately 
\rR/2?ivf for both spins. The larger value of Xr leads to 
two very different wave numbers inside the scattering re- 
gion, fc±, different group velocities for the charge carriers, 
and different indices n±. In the case of XrR/Tivf = 30, 
Figure [3j we notice a clear doubling of the cusps, es- 
pecially evident in panel c for the total wave function. 
Moreover, the cusp positions are in agreement with the 
anticipated values shown by Eq. [I] as one can easily ver- 
ify- 

Let us comment on the system conditions for the ob- 
servation of these effects. We have first assumed a semi- 
classical limit, with kR ^> 1 and k±R ^> 1; second, we 
have assumed the absence of intervalley scattering (K 
to K'\ which requires that the characteristic length of 
potential change must be much larger than the lattice 
constant. These two rather sensible requirements intro- 
duce constraints on the values of the Fermi energy of 
the incoming electron flux, the gate voltage, and the 
size of the gated region. Recent experiments have re- 
ported values of the Rashba interaction in the range 10- 
lOOmeV for graphene samples grown on Ni and inter- 



calated with Au atoms [21]. This coupling is 2-3 or- 
ders larger than the spatially random spin orbit gen- 
erated by intrinsic ripples in graphene [20] . Assum- 
ing possible values of these parameters in experiments 
as E = 80meV, corresponding to carrier concentration 
n c ~ 0.5 x 10 12 cra~ 2 , V = 200raeV, and X R = lOraeV, 
would require R > lOOOnm and lead to n_ — n + > 0.24, 
allowing the observation of the patterns shown in Fig- 
ure [3J3, and its clear detection in scanning probe exper- 
iments, as \x^ usp — x~ usp \ > 42nm. We should also no- 
tice that although we have used spin-polarized injection 
in our calculations, to expose the spin character of the 
caustic pattern, this is not a necessary ingredient to ob- 
serve the birefringence. Unpolarized electron injection 
will obviously erase the net spin structures we have dis- 
cussed. However, the appearance of duplicate caustics 
and cusps would persist and result in an identical spatial 
structure of the total wave function to the ones shown 
above. Similarly, K-K' coupling would preserve the du- 
plicate caustic pattern, although with decreased contrast, 
as the enhanced backscattering would reduce the overall 
transparency of the region. 

In conclusion, we studied the effects of Rashba SOI on 
the scattering of electrons in graphene. For spin polar- 
ized injection, we find a spatial modulation of the carrier 
spin, with characteristics that depend on the strength of 
the Rashba interaction and may be important in the im- 
plementation of spin based devices in graphene systems. 
More important, we find the selective formation of spinful 
cusps and caustics in the region containing large Rashba 
SOI, opening the possibility of spin beam filters and split- 
ters. Moreover, the manifestation of electronic birefrin- 
gence at relatively large but experimentally achievable 
values of the Rashba SOI [21 , is the doubling of caustics 
and cups produced by these refringent Veselago lenses; 
the spacing between the two different chiral cusps is pro- 
portional to the strength of the Rashba interaction in 
the system. As such, the birefringence can be used to 
determine the strength of the Rashba SOI in the region 
of interest in the system. 

We believe that the concept of birefringent electron op- 
tics established here can lead to important consequences 
in experiments and applications. For example, using 
guiding gates [7] over graphene samples grown on Ni[lll], 
or creating a suitable spatial pattern of gold intercala- 
tion, one could design a birefringent waveguide with two 
different critical angles (6± = sin -1 (k±/k)) for the two 
orthogonal chiral states of the system, leading to the 
independent propagation of chiral states, with the cor- 
responding spatial modulation of the spin componenets 
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